Tight-binding study of interface states in semiconductor heterojunctions 
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Localized interface states in abrupt semiconductor heterojunctions are studied within a tight- 
binding model. The intention is to provide a microscopic foundation for the results of similar studies 
which were based on the two-band model within the envelope function approximation. In a two- 
dimensional description, the tight-binding Hamiltonian is constructed such that the Dirac-like bulk 
spectrum of the two-band model is recovered in the continuum limit. Localized states in heterojunc- 
tions are shown to occur under conditions equivalent to those of the two-band model. In particular, 
shallow interface states are identified in non-inverted junctions with intersecting bulk dispersion 
curves. As a specific example, the GaSb-AlSb hetero junction is considered. The matching condi- 
tions of the envelope function approximation are analyzed within the tight-binding description. 



PACS number(s): 73.20.-r, 73.40.Lq, 73.20.Fz 



I. INTRODUCTION 



The advent of highly developed methods of nanotechnology has opened rich possibilities to study quantum effects in 
i-^J \ semiconductor heterostructures, such as heterojunctions, quantum wells, and superlattices (for a recent survey, see Ref. 

jjj). Continued interest is placed on the occurrence and properties of localized electronic states in heterostructures, 
i.e., states that are bound in the growth direction of the structure. This interest arises mainly from the effect that 
i^i - this kind of state is expected to have on transport properties of nanodevices |@||. Basically, two different types of 
localized states can be distinguished. In quantum wells and superlattices, such states occur owing to the presence of 
minima in the band edge potential that confine the electronic motion in the growth direction ("size quantization"). 
On the other hand, the existence of localized states at the interface of a single heterojunction, when no confining 
potential is present, hinges on the coupling between the electronic motion perpendicular to the interface plane and 
y—{ \ the free motion parallel to this plane . Owing to this coupling, energy can be exchanged between perpendicular 
■ and parallel degrees of freedom. 

While localized states in confining potentials have been widely studied [|]]|](|, only few theoretical investigations 
have dealt with localized states at single heterojunction interfaces j|,0-[ll]]. These were carried out within the frame of 
the two-band model [J^|,[l| in the envelope function approximation, using an effective Dirac-type Hamiltonian where 
the interaction between the motion perpendicular and parallel to the interface is mediated by an effective spin-orbit 
coupling. It was found that in abrupt junctions with band inversion localized interface states occur with a linear 
"in-plane" dispersion . In a more general treatment 8 , interface states with a nonlinear in-plane dispersion were 
identified in junctions without band inversion, provided the ratio of energy gap and effective mass changes across the 
interface and the bulk dispersion curves of the two semiconductors intersect. Recently, the two-band model has been 
reformulated in terms of a one-band (Schrodinger- like) description ||. The qualitative features for abrupt junctions 
were found to persist essentially for graded junctions flMPHTlll . 

The two-band model in conjunction with the envelope function approximation allows a simple and transparent 
description of localized interface states. However, a more fundamental approach would have to take explicitly into 



account the atomistic structure of the semiconductors making up the heterojunction. Here, the tight-binding ap- 
proximation, which has been widely used for the description of bulk matter H,[j~5| and surface states Jig], offers 



an appropriate starting point. In the present paper, we apply the tight-binding method to study interface states 
in heterojunctions, with the aim to provide a microscopic foundation for the results obtained from the two-band 
treatment. We formulate a tight-binding model for heterojunctions that is based on a simplified description of III-V 
bulk semiconductors. A two-dimensional Hamiltonian with nearest- neighbor interactions is set up from which the 
Dirac-like bulk spectrum of the two-band model is recovered in the continuum limit and which allows an analytic 
derivation of the conditions for the existence of interface states in heterojunctions. 

The paper is organized as follows. In Sec. II, the tight-binding model for heterojunctions is formulated and a 
relation for the in-plane energy dispersion of localized interface states is derived. The correspondence between the 
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parameters of the two-band model and of the tight-binding model is revealed. In Sec. Ill, the properties of the energy 
spectrum of interface states are analyzed and compared to those of the two-band model; the case of the GaSb-AlSb 
heterojunctions is discussed as a specific example. The matching conditions for the envelope functions of the two-band 
model are elucidated in the tight-binding model. Section IV contains a summary and our conclusions. 



II. TIGHT-BINDING MODEL FOR HETEROJUNCTIONS 



A. Qualitative picture 

We consider heterojunctions composed of two different materials of the III-V class, thereby covering most of the 
cases of practical interest ]l7j . A reasonable description of the band structure of homogeneous III-V semiconductors 
can be achieved within the tight-binding approach , assuming a minimal basis set of one s-orbital and three p- 
orbitals on each atom and including only nearest-neighbor interactions. As we are looking here for localized states 
lying in the fundamental energy gap of the heterojunction (i.e., for states lying between the maximum of the valence 
band edge profile and the minimum of the conduction band edge profile of the junction), we are interested only in an 
accurate description of the bands of the homogeneous semiconductors closest to their respective gap. Then, following 
Ref. pl| where a tight-binding calculation of the band structure of the InAs-GaSb superlattice was performed, we can 
reduce the size of the basis set by discarding the p-orbitals on the metallic atoms and the s-orbitals on the nonmetallic 
atoms. 
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FIG. 1. Two-dimensional square lattice representing a heterojunction formed by two III-V semiconductors. Integers attached 
to the a;-axis and z-axis correspond to labels n and m, respectively. Boxes with dashed lines enclose unit cells with n, m = 0, — 1 
and n,m — 0,0, respectively [cf. Eq. (p^)]; the box with full lines encloses atomic sites entering Eq. (pcj). Symbols in unit cells 
designate s-orbitals corresponding to coefficients u nm (open boxes) and u„ m (shaded boxes), and p-orbitals corresponding to 
coefficients v„m (open circles) and v„ m (shaded circles) [cf. Eq. (^|)]. On-site energies are eF and e^ ] ; j = 1 for m < (z < 0), 
j = 2 for m > (z > 0), and the hopping matrix elements between nearest neighbors are ±Vj. Hopping across the interface at 
z — is described by matrix elements ±V. 



For the actual treatment of the heterojunction, we consider a simplified, two-dimensional model system that allows 
a transparent description in which the properties of interface states can be derived essentially in analytic form. The 
crucial feature of this model system is the explicit inclusion of the coupling between the electronic motion perpendicular 
and parallel to the interface plane. Basically, a two-dimensional description appears adequate for the present problem 
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as in the continuum limit one deals, due to axial symmetry about the direction perpendicular to the interface, with 
a truly two-dimensional problem. 

As depicted schematically in Fig. [I], we consider an (infinitely extended) square lattice with lattice constant a, where 
the z-axis is chosen perpendicular to the interface plane (with the interface at z = 0) and the x-axis parallel to it. 
With s- and p-orbitals assigned alternatingly to the atomic sites, the unit cells of the model system are each formed 
by two s-atoms and two p-atoms. Integers n and m corresponding to the x-direction and the z-direction, respectively, 
are used to label the position of the unit cells on the lattice. The homogeneous semiconducting materials 1 and 2 
forming the hetero junction are characterized by atomic energies ei^ and e p ^ (j — 1,2), and by the magnitude Vj of 
the hopping matrix elements between nearest-neighbor atoms, which is assumed to be independent of the hopping 
direction. As the p-orbitals are odd, the matrix elements connecting a specific atomic orbital to its nearest neighbors 
in either of the two directions differ in sign. The interface separating the semiconductors is located between the 
vertical rows of units cells with labels m = — 1 and m — 0. A value V is assumed for the magnitude of the matrix 
elements connecting s- and p-orbitals across the interface; in the course of the formulation of the model, V will be 
expressed in terms of the bulk matrix elements V\ and Vi . 

In the model system illustrated by Fig. [j], the relative signs of the hopping matrix elements in adjacent linear 
chains of atoms have been chosen such that the matrix elements connecting the atoms within a unit cell add up to 
a nonzero value. By reversing in Fig. [I] the signs of all matrix elements within every second chain either in the x- 
or the z-direction, one would obtain a topologically different system in which the sum of the matrix elements within 
a unit cell is zero. The distinctive feature of the model system of Fig. |l| is that in the continuum limit it exhibits 
a Dirac-like bulk dispersion (cf. below). It is noted that there is a close analogy between the present lattice model 
and certain two-dimensional models describing "flux phases" in strongly correlated electron systems, like high-T c 
superconductors Jl9[-p2|. 



B. Homogeneous semiconductors 

Turning now to the detailed formulation of our model, we write the two-dimensional tight-binding Hamiltonian for 
a homogeneous semiconductor characterized by atomic energies e s and e p and by the magnitude V of the nearest- 
neighbor hopping matrix elements as 



H ^ ] { e s( a lim a nm + a lim a nm) + £ p(^nm^nm + ^nmPnm) 

+ V\(Pn,m— 1 + ^n— l,m — ^nm ~ ^nm) a nm + (^n+l,m ~ ^n,m+l + — Kim^nm 

1 V n,m+l 1 n— l,m ^nm nm) nm 1 V^n+l,?™ n,m — 1 1 nm ^nmJ^nml } \ J 

(cf. Fig. Q). Here, the operators a nm ,a'l im (a nm ,a nm ) and b nm ,b nm (b nm ,b nm ) create (annihilate) electrons in s- and 
p-orbitals, respectively, in the unit cell labelled n, m. The energies e s and e p do not depend on n and m. Representing 
the one-electron wave function $ as a linear combination of s- and p-orbitals, 

^ = ^ ^ { u nm a nm F ^nm ^tim F v nm b\ im + V n m b nrn )\0) ' (2) 

n,m 



we write the Schrodinger equation {H — E)^ = in terms of a coupled set of equations for the coefficients u nm , u nm . 
Vnrm and v nm as 





= o, 




= o, 




1=0, 




1=0. 



(e s - E)u nm + V(v ntm -x + v n -i, m - v nm - 
(e s - E)u nm + V(v n+ i tm 

(^p -^)^nm ~t~ ^(^n.m+1 ~t~ ^n-l,m ^nm 

( e P _ E)v nm + V(u n+ i tm - U n ,m-1 + U nm ~ U n m) = . (3) 

Exploiting the periodicity of the model system, we can find a general solution of Eqs. (||) in the form of Bloch sums 
with coefficients 

u nm = uqq enp[i(n(f) x + m(f> z )} , (4) 

and similarly for u nm , v nrn , and v nm . The phases <j) x and <f> z are related to the wave numbers k x and k z through 
4> x = k x a and (f> z — k z a. Insertion of Eqs. (||) into Eqs. (||) leads to the system of equations 
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with Z = 1 



b ' and X = 1 



been taken at midgap, i.e., e s 



(A - E)uoo - V(Zv o +Xv o ) = 0, 
(A - £)u 00 + V(Z*tfoo -^*w o ) = , 
(-A - E)v 0Q - V(Z*u o+Xu 00 ) = , 
(- A - E)v 00 + V(Zu 0Q -X*u 0Q ) = , 



(5) 



. The quantity A = (e s — e p )/2 is half the energy gap, and the energy zero has 
A (cf. left half of Fig. @). 
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FIG. 2. Schematic band edge profile for an abrupt heterojunction. 



From Eqs. (|B|), the energy dispersion of the homogeneous system is obtained as 

E 2 = A 2 + W 2 [sm 2 ((j} x /2) + sin 2 (<^/2)] . 
A set of four linearly independent (unnormalized) eigenvectors of Eqs. (|^) is given by 



MOO 



1 , woo = , woo = - 



VZ* 



and 



"oo = , uqo = 1 , voo = 



A + E ± 



VX 



E± 



woo 



woo 



VX* 
'A + E± 



VZ 



E± 



(6) 



(7) 



(8) 



where E+ (E-) is the positive (negative) root of Eq. (g). 

Equation (^|) describes the dispersion of two adjacent bands in the two-dimensional tight-binding model for homo- 
geneous semiconductors: the positive and negative root correspond to the conduction band and the valence band, 
respectively. In the continuum (long-wavelength) limit k x a, k z a — > 0, we have 



E 1 



A 2 + V 2 a 2 (fc 2 + fc 2 ) 



(9) 



The Dirac-like bulk spectrum of the two-band model is now recovered if the hopping velocity w = Va/h is identified 
with the interband velocity matrix element v [Q . Similarly, the set of eigenvectors ( |7|) a nd (|J) reduces in the continuum 
limit to the usual set of spinor solutions of the free Dirac equation (cf., e.g., Ref. ]23|). The tight-binding model for 
homogeneous systems (in which relativistic effects, in particular spin degrees of freedom, are not explicitly included) 
thus turns out to be equivalent, in the continuum limit, to the relativistic free-electron theory. 

It is worth mentioning that if the signs of the hopping matrix elements were chosen such that their sum within a 
unit cell is zero (cf . above) , one would obtain, instead of Eq. (^|) , a dispersion of the form 



E = A 



W 2 [sin(^/2) ±sin(,/» z /2)] 5 



(10) 



In the present context, this dispersion is ruled out owing to its inappropriate behavior in the continuum limit. Whether 
lattice models with dispersions of the form (10) are useful for other applications appears to be an open question. 
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C. One-dimensional heterojunction 



Having established our two-dimensional tight-binding description of homogeneous semiconductors, we now consider 
heterojunctions and derive the conditions for the existence of localized interface states. Some general features of this 
problem can be elucidated already by looking at the one- dimensional system represented by a linear chain of atoms 
oriented perpendicular to the interface. 

For definiteness, we single out the chain formed by the bottom atoms in all unit cells with n = in Fig. [l]. The 
unit cells of this chain each contain two atoms and are labelled by the index m (we omit the index n = 0). The Bloch 
coefficients for the homogeneous one-dimensional system have the form 



u exp(im4> z ) , v m = v exp(zra</> 2 ) 



(11) 



Localized wave functions falling off exponentially on either side of the heterojunction are characterized by purely 
imaginary phases in the Bloch coefficients corresponding to the homogeneous semiconductors 1 and 2. Writing 



4 13 = 



-lOi 



>o) , 



we obtain from Eqs. ([n]) the relations 

U-i exp[(m + 
Mo exp(— m6 2 ) 



u, 
u, 



1)01] 



u_i exp[(m + 
vq exp(-rn# 2 ) 



for m < : 
for m > 



(12) 



(13) 



In order to connect the decaying wave functions across the interface, we write down the set of coupled equations for 
the coefficients uq, Vq corresponding to the atomic orbitals in the unit cells labelled m = — 1 and m = that 

sandwich the interface. Using the relations ( p"3| ) to express the coefficients v^ 2 and u\ in terms of the coefficients u_i 
and uq , respectively, we have from the first and third of Eqs. (|^) 



(Ai - E)u-i - V\Z x v- X 

("Al-B)!)-! -VlU_! 

(A 2 + W - E)u - V 2 v 

(-A a + W- E)v - V 2 Z 2 u 



+ V'u Q 



(14) 



where 



1,2 



The quantities Aj — (ei^ — ep^)/2 are half the energy gaps, and W — (ef + e^')/2 is the work function offset (the 
energy zero is taken at midgap in semiconductor 1; cf. Fig. ||). 

Eliminating in the equations involving the matrix element V' the coefficients u—± and vq by means of the first and 
fourth equation and using the bulk energy dispersions [cf. Eq. (^) specialized to the one-dimensional case] for the 
imaginary phases given by Eqs. (p"2|), we obtain from Eqs. (Jl4|) a system of two equations for uo and V—i. For its 
solution, we must require 



(15) 



V?Vi(e ei - l)(e 92 - 1) + U' 2 (A! - E)(A 2 - W + E) = 



(16) 



It is seen from Eq. (H) that for any choice of the values of the parameters V\, V%, and V', there are no solutions 
within the energy range max{— Ai, — (A 2 — W)} < E < min{Ai, A 2 + W}, i.e., within the fundamental energy gap 
of the heterojunction. We thus find that there are no localized (bound) interface states in the one-dimensional case, 
a result that was to be expected since here no energy can be stored in the motion parallel to the interface. 



D. Two-dimensional heterojunction 

We now turn to the two-dimensional tight-binding treatment of heterojunctions. Here, localized states are described 
by two-dimensional Bloch waves with real phase (j) x — k x a corresponding to free motion in the direction parallel to 
the interface and purely imaginary phases <pz given by Eqs. (|l2] ) which correspond to exponential fall-off in the 
z-direction on either side of the interface. From Eqs. (Q), we then have 
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"o,-i exp[iruf> x + (m 4- l)#i] 



(17) 



for m < 0, 



"wn = "oo exp[m<^ x - m6 2 



(18) 



for m > 0, and similarly for u nm , v nm , and v nm . 

As in the one-dimensional case, we now write down the set of coupled equations for the coefficients corresponding 
to the atomic sites in the unit cells labelled n, m = 0, — 1 and n, m = 0, (boxes with dashed lines in Fig. |l]). Using 
relations ( P~7[ ) and (18) to eliminate all coefficients with labels other than n, m = 0, —1 or n,m = 0,0, we have from 
Eqs. © 



(Ai -£> ,_i 
(-A 1 -E)v 0> _ 1 
(Ax --E)uo,_i 

(-A x - -E)«0,-i 
(A 2 + - £>oo 
(-A a + W - 
(A 2 + IU - £)w 00 
(-A a + W - £> 00 



Tq(Ziu ,-i 
Vi(^iUo,_i 
Vi(uo,-i 
Vi(uo,-i 
■ Va(«oo 

^("00 

V 2 (Z 2 v 00 

V 2 (Z 2 UQ0 



+ Xv ,-i 

- **Uo,-l 

- **«0,-l 

+ Im 0) -i 
+ -^"oo 

- X*u 00 

- X*v 00 
+ Xuqo 



V'v 00 
V'u oa 
V'v ,-i 
V'u ,-i 



(19) 



The coefficients "o,-i, vo,—i, "oo; and i>oo in the four equations involving the coupling V' can be eliminated with the 
help of the remaining equations, leading to a closed system of equations for the coefficients Wo,-i, "o.-i; u oo, and "oo 
corresponding to the four atomic sites immediately adjacent to the interface (box with full lines in Fig. ^): 



(Ai -E)uq,-i 
{-A 1 -E)v 0i -i 
(A 2 + W - E)u 00 
(-A 2 + W - £)v 00 



ViX* Uo ,_i - V'Z^oa = , 

ViXUa,-i + V'Zmoo = , 

V 2 Xv 00 + V'Z 2 v ,-i = , 

V 2 X*u 00 - V%fio,-i = . 



(20) 



Setting the determinant of this system equal to zero and making use of the bulk dispersions of the semiconductors 1 
and 2, 



E 2 =A\ 
{E-Wf = A 2 

[cf. Eqs. (|) and (Q], we obtain the relation 



V(V 2 \e^ ~\)(e^ -l) + V'\l 



4U 1 2 [sin 2 (0 x /2)-sinh 2 (0 1 /2)] 
4U 2 2 [sin 2 (0j2)-sinh 2 (0 2 /2)] 



) + 2y' 2 {AiA 2 +AViV 2 sm z {(j) x /2) - E{E - W)} = 



(21) 



(22) 



Here, the quantities 0j can be expressed in terms of the energy E and the phase <j) x = k x a through Eqs. (21). Relation 
( p2| ) thus determines in an implicit way the in-plane dispersion E(k x ) for localized interface states in a semiconductor 
heterojunction. It represents a principal result of the present paper. 

Up to this point, the matrix element V' describing hopping across the interface has been considered an independent 
parameter. In order to fix the value of V within the present model in a consistent way, we relate it to the bulk matrix 



elements V\ and V 2 by considering the behavior of the wave function coefficients u 
interface is crossed. 



and v nm when the 



For u nm , e.g., this behavior is readily inferred by combining the first and second of Eqs. (JlS|) , the second of Eqs. 
(p0|), and the first of the bulk dispersions (21). With similar reasoning for the other coefficients, we find 



and 



"00 
"0,-1 



"00 

"0,-1 



"00 



= V Le e 
"o,-i V 



"00 
"0,-1 



YL 
v 2 



(23) 



(24) 



Interpreting the set of coefficients referring to a specific unit cell as the components of a single, four-component wave 
function, we require all coefficients to change across the interface by one and the same factor, i.e., we equate the 



right-hand sides of Eqs. (|23|) and (p4T). This determines U' 2 as 



G 



V' 2 = VxV 2 e e ^ , 



(25) 



and all coefficients change by the factor (V1/V2) 1 ' 2 e^ 6l ~ S2 '' 2 . In the continuum limit 6\, 6 2 — > 0, the matrix element 
V' thus appears as the geometric mean of V% and V 2 . It is seen below that this choice for V' is indeed necessary in 
order to recover the two-band model. 



III. ANALYSIS AND COMPARISON WITH TWO-BAND MODEL 



A. Energy spectrum 

1. General properties 



Equation ( p2[ ) in conjunction with Eqs. (21) reduces to a quadratic equation for sin (<p x /2) as a function of the 
energy E from whose solution one obtains the in-plane dispersion E(k x ). Mere inspection of Eq. ( p2| ) already provides 
us with conditions for the existence of solutions for 6x, 2 > 0, i.e., conditions for the existence of localized interface 
states. 

It is immediately seen that the spectrum E{k x ) is bounded by the condition 

E(k x ) > E(k x ) , (26) 

where E(k x ) is defined by 

E(E -W) = AxA 2 + WxVi sin 2 (^/2) (27) 

[we consider here only the positive (conduction band) branches of the dispersions for systems without band inversion, 
i.e., systems with Ai, A2 > 0]. The function E(k x ) can be viewed as a kind of mean in-plane bulk dispersion with 
energy gap and hopping matrix element given by the geometric mean of the gaps and matrix elements, respectively, 
of the bulk semiconductors 1 and 2, whose in-plane dispersions £1,2 (Mr) follow from Eqs. (|2|) by setting 9\ = 0% = 0: 

E\ = A 2 + 4V?sin 2 (<fe/2) , 
(E 2 - Wf = A 2 . + W 2 2 sin 2 (<^/2) . (28) 

Moreover, from Eqs. (pl| ) with 9i,9 2 > 0, we have 

E{k x ) < E 1 (k x ) , E{k x ) < E 2 {k x ) . (29) 

For localized states to occur, it is hence necessary that E(k x ) lies below the bulk dispersions of both semiconductors. 
On the other hand, condition ( ^6| ) requires E(k x ) to lie above the mean dispersion E(k x ). Conditions ( ^6|) and ( p9|) can 
be simultaneously fulfilled only if the bulk dispersions E\(k x ) and E 2 (k x ) intersect. The latter condition is equivalent 
to that of the two-band model . 

From Eqs. (E^), it is seen that in the special case where W — the bulk dispersions intersect if and only if the 
condition 

< < 1 (30) 

" AV+V- ~ v ' 

holds, where A± — Ai ± A2, V± = V\ ± V 2 . This implies, in particular, that A_/V^_ < 0. If the bulk dispersions 
intersect, the mean dispersion lies between the bulk dispersions for small values of k x (cf. Fig. |^ below), but below 
both bulk dispersions in the vicinity of the point of intersection (thus the attribute "mean" given to the function E is 
not to be interpreted too narrowly). At k x — n/a, where the bulk curves and the mean curve all attain a maximum, 
the mean curve lies between the bulk curves if the conditions 



hold simultaneously for j — 1 and j = 2. In this case, the energy spectrum of localized interface states, E{k x ), is 
restricted to a finite fc^-range centered around the point of intersection of the bulk dispersions, with end points kx 
(j = 1, 2) determined by the points of intersection of E(k x ) with E\{k x ) and E 2 (k x ). For W = 0, we have 
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4 J ' } = - arcsin(-A j A_/[4V r jF-]) 1/2 



(32) 



Note that the end points of the spectrum do not depend on the matrix element V', i.e., on properties characterizing 
the interface, although the dispersion determined by Eq. ( p2|) depends on V . 

With V' expressed through Eq. (p5j), Eq. ( p2j ) can be easily solved in the continuum limit, keeping terms of lowest 
order in 9\ , 62 and <j) x . The resulting energy spectrum is 

^(M = ,2,, 9 7T {^(?i 2 fc> lW _ + AiA_) ± k x d{h 2 k 2 x w 2 _ + A 2 _ - W 2 ) 1 ' 2 } , (33) 

where Wj = Vja/h are the hopping velocities, u>_ = w\ — W2, and d = W1A2 — u^Ai. This spectrum agrees with that 
of the two-band model (Eq. (5) of Ref. Q ) if the identification 

Wj ee vj (34) 

is made, with Vj the interband velocity matrix elements of the two-band model. 

The tight-binding model thus generalizes the two-band model and provides a microscopic foundation for the latter. 
In particular, the shallow interface states identified in Ref. [Q for the case of non-inverted semiconductors with 
intersecting bulk dispersion curves thereby find a firm theoretical basis. 



2. Example: GaSb-AlSb 

In order to apply our tight-binding model to specific heterojunctions, we have to fix the values of the parameters Ai, 
A2, W, Vi, and V%. For junctions formed of common III-V materials, the values of the first three of these parameters 
can be taken from experiment (cf., e.g., Ref. [jl^,^4|). For the hopping matrix elements Vj, different sets of values are 
available |p^ , ^5| , ^5|j26| cither from theoretical estimates or from fits to specific material properties. It is not obvious 
which one of these sets represents an adequate choice in the context of the present, simplified model. We therefore 
resort here to expressing Vj in terms of the velocity matrix elements Vj of the two-band model as Vj = hvj/a (cf. Sec. 
III.A.l). The velocity matrix elements, in turn, are given Q by Vj — (Aj /m*) 1 / 2 , where m* is the conduction-band 
effective mass at the T-point. 

As a specific example, we consider the conduction band states in the GaSb-AlSb heterojunction (cf. Fig. |J). With 
the values for Ai, A 2 , and W taken from Ref. Q and the (theoretical) values for m\ and taken from Ref. [ jl5| , 
we find that the in-plane bulk dispersion curves E\ and E2 corresponding to GaSb and AlSb, respectively, intersect 

at <p x — 1.526. The mean in-plane bulk dispersion curve E intersects the curves E\ and E% at (f>^ = 1.200 and 
(2) 

4> x = 1.874, respectively. Localized interface states exist within the ^-range between these limits, with the in-plane 
energy spectrum shown in the inset of Fig. ^. The binding energy of these states (as determined by the difference in 
energy to the bulk dispersion curves) is of the order of 10 meV in the vicinity of the point of intersection of the bulk 
curves. Since the Dirac-like dispersion curves of the two-band model intersect at a point close to the intersection of the 
tight-binding curves (cf. Fig. |^), the results of the two-band model do not differ much from those of the tight-binding 
model. 

From considering the example of Fig. ^|, it becomes apparent that the properties of the calculated energy spectrum 
of interface states depend sensitively on the values of the model parameters Aj , W, Vj . Fairly small changes in one of 
these values may lead to a large shift in the position of the point of intersection of the bulk dispersion curves, or may 
even remove the intersection. With regard to the experimental verification of interface states, it would, therefore, be 
premature to attach quantitative significance to the present results. In order to achieve a fully reliable description of 
interface states, refined tight-binding calculations including the determination of a consistent set of parameters (as 
outlined, for example, in Refs. |p7| , p8[ ) are called for. 



B. Wave functions 

In the tight-binding approach, the wave function ^ represented by the expansion (g) and determined as a solution 
of Eqs. ( p0| ) describes the electronic motion in the heterojunction as a whole. The change of the wave function across 
the interface is uniquely fixed by the change in the coefficients expressed by the relations ( |23| ) and ( p4| ) with V' given 
by Eq. (§]). 
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FIG. 3. Dispersion curves for the conduction band states in the GaSb-AlSb heterojunction, calculated with the parameter 
values Ai = 0.365 eV, A 2 = 0.79 eV, W = 0.075 eV (taken from Table 2 of Re£_M), Vi = 1.3617 eV, V 2 = 1.2041 eV 
(derived from the effective masses given in the second column of Table 6.2 of Ref. [p_5j; cf. text). Solid curves: in- plane bulk 
dispersions E\, Ei given by Eqs. Short-dashed curve: mean in-plane bulk dispersion E given by Eq. (|2^). Long-dashed 

curves: Dirac-like in-plane bulk dispersions of the two-band model [cf. Eq. (^)]. Thick solid curve in the inset: in-plane energy 
spectrum of localized interface states, obtained from the numerical solution of the tight-binding equations ( [ill ) and ( p^ ) with 
V' expressed through Eq. (^H|). The arrows in the inset mark the intersection of the bulk dispersions and the intersections of 
the mean dispersion curve with the bulk dispersions, respectively. 



In the two-band model for abrupt heterojunctions, on the other hand, separate solutions are obtained for the wave 
functions in the two bulk regions of the junction, which are subsequently matched at the interface. The role of the 
matching conditions has been analyzed in Ref. Q (cf. also Ref. p9| ), where conditions for the physical envelope 
function were established which relate the discontinuity of the wave function to that of the interband velocity matrix 
element across the interface. 

In the following, we demonstrate the consistency of the two-band and tight-binding descriptions with regard to 
the behavior of the wave functions across the interface. We consider the continuum limit of the system of equations 
(|J) describing the homogeneous semiconductor; omitting the labels attached to the coefficients Moo > ^oo , ^oo , ^oo , we 
introduce linear combinations 

ip\ = u + iXu ; xx — v + i (35) 

where A = ±1. In the continuum limit, keeping terms up to first order in cf> z and 4> x and considering ip\ and x\ to be 
continuous functions of z and x, i.e., ip\(z, x), X\{ z i x) oc exp[i(k z z + k x x)\, we replace the wave numbers k z and k x 
with — id z and —id x , respectively, so that Eqs. @ can be expressed in the form 

(A - E x )^ x {z, x) - Va(d z - i\d x ) X x(z, x)=0, 

(36) 

Va{d z + i\d x )ip x (z, au) - (A + E x ) X x(z 1 x)=0. 

These equations are equivalent to the Dirac-like bulk equations of the two-band model (Eqs. (2) of Ref. Q) if the 
hopping velocity w — Va/fi is replaced with the interband velocity matrix element v [cf. Eq. (^J)] and if the parameter 
A is identified with the eigenvalue of the helicity operator. 

Assuming the interface to be located at z = and using Eqs. fl23|), (p4|), and (p5|), we can now relate the values 
= ip\(z = —0) and x x = X\{ z — — 0) to the values tp^ =ip\{z = +0) and X\= X\( z = +0) : 



;(!) (!) /tr \ V2 



^ xi 2) \Vi 



(37) 
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These relations are seen to be equivalent to the matching conditions of the two-band model (Eqs. (3) of Ref. if 
we adopt Eq. ( |34| ) to relate the hopping matrix elements Vj to the interband velocity matrix elements Vj . 

IV. SUMMARY AND CONCLUSIONS 

We have presented a tight-binding study of localized interface states in semiconductor hetero junctions. Starting 
from the usual description of homogeneous III-V semiconductors, we have set up a two-dimensional lattice model, in 
which the coupling between the electronic motion perpendicular and parallel to the interface plane in a junction is 
explicitly taken into account. Yet, this model is simple enough in order to allow an essentially analytic derivation of 
the conditions for the existence of interface states and a transparent discussion of their properties. 

In constructing the model Hamiltonian, an ambiguity arises in the choice of the relative sign of the hopping matrix 
elements in adjacent linear chains of atoms. Depending on this choice, two topologically different sytems are obtained 
in which the matrix elements within a unit cell add up either to zero or to a nonzero value. In the present context, 
we choose the alternative with nonzero sum, which leads in the continuum limit to the Dirac-like description of 
homogeneous semiconductors within the two-band model. 

As a principal result, we have derived an implicit expression for the in-plane energy dispersion of heterojunctions. 
From this expression, conditions for the existence of localized interface states are inferred, which are equivalent to 
those of the two-band model. In the continuum limit, the energy spectrum of the two-band model is recovered if the 
hopping velocities of the tight-binding model are identified with the interband velocity matrix elements of the two- 
band model. By analyzing the behavior of the tight-binding wave functions across the interface, we have elucidated 
the matching conditions for the physical envelope function in the two-band model. 

In conclusion, our tight-binding model for interface states in semiconductor heterojunctions generalizes the two-band 
model and provides a microscopic foundation for it. Although our model is based on the tight- binding description 
of III-V semiconductors, its qualitative features do not appear too sensitive to specific details of the model and its 
results are therefore assumed to be valid for broader classes of materials. Considering the possibility of an experimen- 
tal verification of the type of interface state identifed in the present work, we have emphasized the need for refined 
tight-binding calculations which go beyond the frame of our simple model and in which consistent parameter sets are 
used. 



Note added in proof: After acceptance of this manuscript, Ref. ]3C| ] was called to our attention. In that paper, 
interface states in heterojunctions are treated in a one- dimensional tight-binding approach that uses, with respect to 
the description of III-V semiconductors, more general assumptions than our two-dimensional model. A relation for 
the energy of interface states is obtained that generalizes our relation (|l6| ) and for specific parameter values admits 
solutions within the fundamental energy gap of the junction. As the electronic motion parallel to the interface is 
disregarded in Ref. |30|, that paper does not deal with the in-plane dispersion which is the subject of the present 
work. 
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